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Abstract 

In this paper, we consider a model in which the unlicensed or the Secondary User (SU) equipped 
with a Cognitive Radio (CR) (together referred to as CR) interweaves its transmission with that of 
the licensed or the Primary User (PU). In this model, when the CR detects the PU to be (i) busy 
it does not transmit and; (ii) PU to be idle it transmits. Two situations based on CR's detection of 
PU are considered, where the CR detects PU (i) perfectly - referred to as the "ideal case" and; (ii) 
imperfectly - referred to as "non ideal case". For both the cases we bring out the rate region, sum 
capacity of PU and CR and spectral efficiency factor - the ratio of sum capacity of PU and CR to 
the capacity of PU without CR. We consider the Rayleigh fading channel to provide insight to our 
results. For the ideal case we study the effect of PU occupancy on spectral efficiency factor. For the 
non ideal case, in addition to the effect of occupancy, we study the effect of false alarm and missed 
detection on the rate region and spectral efficiency factor. We characterize the set of values of false 
alarm and missed detection probabilities for which the system benefits, in the form of admissible 
regions. We show that false alarm has a more profound effect on the spectral efficiency factor than 
missed detection. We also show that when PU occupancy is small, the effects of both false alarm 
and missed detection decrease. Finally, for the standard detection techniques viz. energy detection, 
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matched filter and magnitude squared coherence, we show that that the matched filter performs best 
followed by magnitude squared coherence followed by energy detection with respect to spectral 
efficiency factor. 

I. Introduction 

While certain parts of the frequency spectrum are crowded by users, most part of the 
spectrum still remains heavily unoccupied [lj. Due to this spectrum imbalance, the legacy 
command and control policy of the FCC and the recent proliferation in usage of wireless 
devices, new techniques that can perform Dynamic Spectrum Access (DSA) are needed. The 
so called Cognitive Radio (CR), in essence, performs DSA by sensing the radio environment, 
choosing the best channel (frequency, time slot or codes) that is available and adapting 
parameters to communicate simultaneously ensuring the interference to legacy users below a 
prescribed level [|2). 

In literature (for example [|3), flU) there are three paradigms for CR based systems men- 
tioned, based upon the nature of operation of the CR, viz. (i) interference avoiding or 
interweave, (ii) interference limited or underlay and (iii) interference mitigating or overlay. 
Each of the above paradigm assumes different amounts or levels of side information that the 
SU has about PU transmission. 

From an information theoretic viewpoint, most of the recent works focus on modelling the 
CR channel as an interference channel with some side information about the PU codebooks, 
messages, channel characteristics etc, being available to the CR transmitter @)-[15]. Also 
they mostly address underlay and overlay based CRNs. 

In 121 a two switch model for CR is presented with one switch at the SU transmitter and 
the other at the SU receiver. The role of each switch is to identify transmitting opportunity 
(channel) which can be different at the SUs transmitter and receiver. Transmission is then 
done on a channel which is identified to be free on both the transmitter and the receiver ends. 
In H the CR channel is modeled as two sender two receiver model with the transmitters 
having non causal information. In [7],[8] capacity for low and high interference regimes is 
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analyzed, wherein the PU codeword is assumed to be known to the CR. Extending from 
two user case to the three user case, in [|9) achievable rate regions and outer bounds are 
derived. Channels are modeled as interference channels where the transmitters cooperate in a 
unidirectional manner via 3 different noncausal message-sharing mechanism specified therein. 

For the case of underlay based systems, in [fTOl . ffTD the capacity gains offered when only 
partial channel information of the link between the CR transmitter and PU receiver is available 
to the CR under average received power constraint is studied. In |[T2l an optimal power 
allocation policy to maximize capacity, which exploits a two-dimensional frequency-selectivity 
on CRs own channel as well as CR to PU channel under an interference-power constraint 
as well as a conventional transmit power constraint is proposed. Moving further ahead lfT3~l 
proposes a joint overlay and underlay technique to maximize the total capacity of CR. 

While most of the above works either attempt to characterize the achievable regions or 
maximize CR capacity, certain works study about the sum capacity of the CR and PU together 
so that capacity gains of entire system is brought out lfl4l , lfT5l . 

Significantly, all the above works assume that the CR and PU have perfect side information 
about the PU occupancy. Although such models are of use, practical scenarios wherein there 
is imperfect side information available needs to be studied as capacity can decrease with 
imperfect side information, significantly. 

In (H, a simple rate region involving rates of PU and CR has been brought out for the 
case where the CR interweaves its communication with the PU by means of time sharing. It 
is assumed therein that CR has perfect knowledge of PU's channel occupancy. We refer to 
this as the "ideal case". 

In this paper, we find the spectral efficiency factor, which is the ratio of the sum capacity 
of CR and PU to the capacity of the PU in absence of CR, for the ideal case. With the help 
of an example of Rayleigh flat fading channel, we study its variation with PU occupancy for 
different relative power levels of the CR and PU. Along with the ideal case, we consider the 
case where CR interweaves its communication with PU with "only an imperfect" knowledge 
of the PUs channel occupancy. We refer to this as the "non ideal case". We obtain rate regions 
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for the non ideal case and compare it with the ideal case. We show that the rate region of 
the non ideal case is a sub-region of the ideal case. We also discuss the effects of false alarm 
and missed detection on the rate region for the non ideal case. We then find the maximum 
sum capacity and use this expression compute the spectral efficiency factor. We bring out 
effects of false alarm and missed detection on the spectral efficiency factor. We show that, 
probability of false alarm has more profound effect than the missed detection on spectral 
efficiency factor. 

In the ideal case the, spectral efficiency factor is always greater than unity. However, in 
the non ideal case, if the error in detection is large, i.e. the probability of false alarm, pf a , 
and the probability of missed detection, p m d, either individually or collectively transgress 
certain limits; the spectral efficiency factor turns out to be less than unity. We characterize 
regions made up of those (pf a ,Pmd) points, called the admissible regions, for which we obtain 
spectral efficiency factor greater than unity. Finally we compare known detection schemes 
from a spectral efficiency perspective. 

The rest of the paper is as follows. In Section [TT] we discuss the system model. In Section 
ITTll we show the rate region for the case where the CR detects the PU perfectly. We then 
address the non ideal detection and characterize the corresponding rate region in Section [IV] 
We also find conditions for obtaining maximum sum capacity therein. Then we discuss the 
spectral efficiency factor in Section |Vj wherein we discuss the effects of false alarm and 
missed detection on the spectral efficiency factor. We then show a comparison of the certain 
standard detection techniques in light of the results obtained in section |VIJ Finally conclusions 
are given in section IVII1 

II. System Model 

In this paper we assume that the Primary User (PU) transmits with a probability 1 — p and 
the channel is free with a probability p. However, from a detection perspective, the occupancy 
of a single channel is defined as the unconditional probability that the measured signal strength 
does exceeds some predetermined threshold [16|. In spectrum sensing, the detection involves 
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comparing a test statistic against a threshold, suitably chosen. Let Y be the observed random 
variable and X the transmit random variable, so that X = 1 if there is no transmission and 
otherwise. It follows that p = F[X = 1]. We now discuss the system models for the two 
cases based on the relation of Y and X. 

A. System Model for Ideal Detection 

In case of perfect detection, Y = X. This is referred to as the case of "ideal" detection. 
The system model is described in Fig. \T\ 

The baseband relations for the PU and the CR are defined as 

For, Y = X = 0; Y p = h p X p + Z p , 

For, Y = X = 1; Y c = h c X c + Z c , (1) 

where X p and X c are the signals transmitted, Y p and Y c are the signals received by the PU 
and the CR respectively. Z p and Z c are uncorrected circularly symmetric zero mean variance 
a 2 white Gaussian random variables i.i.d. over time for the PU and CR respectively, denoted 
as, Z p , Z c ~ CjV(0, a 2 ). The average power constraints on the PU and CR are ¥.[X P ] < P p 
and E[X C ] < P c . The fading coefficients h p and h c are uncorrelated and i.i.d over time. These 
channels are asymmetrically time shared, i.e. the CR channel comes into existence only when 
the PU channel is not used. 

B. System Model for Non Ideal Detection 

In general, 7^1, which is referred to as the case of "non ideal" detection. In non ideal 
detection there are errors in detection. The probability of false alarm, Pf a — P[V = l\X = 0] 
is the probability that CR detects the channel to be free given it was occupied by the PU. 
Similarly, the probability of missed detection, p md = F[Y = 0\X = 1], is the probability that 
CR detects the channel to be occupied given it was free from PU transmission ffTTl . Note 
that the other two conditional probabilities can be expressed in terms of pf a and p mc [. The 
system model for non ideal case is as shown in Fig. |2l 
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Due to the possibility of false alarm, the PU and CR might transmit at the same time, 
which leads to a possibility of the interference channel. Hence, the channel can be a pure flat 
fading channel with probability ¥[X = 1, Y — 1] and an interference channel with probability 
F[X — 0,Y — 1]. The entire channel is time shared between a pure flat fading channel for 
PU, a pure flat fading channel for the CR and finally an interference channel for the PU and 
CR. We thus have a time sharing random variable as follows. 

Definition 1. Time sharing random variable Q is as defined 

1) Q = 1 is the event that X — 1,Y — 1, CR correctly detects the channel to be free; 

2) Q = 2 is the event that X = 0, Y = 1, CR erroneously detects the channel to be free ; 

3) Q = 3 is the event that X — 1,Y — 0, CR erroneously detects the channel to be 
occupied; 

4) Q = 4 is the event that X = 0, Y = 0, CR correctly detects the channel to be occupied. 
The probability distribution on Q is as follows 



P(Q 


= 1) 


= ¥[X = 


1,Y = 


1] 


= p(l-Pmd), 


P(Q 


= 2) 


= ¥[X = 


0,Y = 


1] 


= iX-p)Pfa, 


P(Q 


= 3) 


= ¥[X = 


1,Y = 


0] 


= PPmd, 


P(Q 


= 4) 


= ¥[X = 


0,Y = 


0] 


= (l-p)(l- P/o ) 



Note that Q = 1,2,4 respectively denote the instances where there is a flat fading channel 
for the CR, an interference channel between the PU and the CR, and a flat fading channel 
for the PU. Also note that Q = 3 is the instance where the channel goes waste as there is 
no transmission during this time, neither from PU nor from CR. The distribution of the time 
sharing random variable is fixed for a given (pf a ,Pmd)- We have the baseband representation 
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as 



For Q 



1: 



Y c — h c X c + Z, 



For Q 



2; 



Y c — h c X c + h pc X p + Z c and Y p — h p X p + h cp X c + Zj 



For Q 



4; 



(3) 



where Yp, Y" c , X p , X c , Z p , Z c , h p and h c remain the same as in the ideal case. h pc and h cp 
are also i.i.d over time and uncorrelated with each other as well as with h p and h c . 

The system model for the non ideal case is a very general model and encompasses some of 
the models available. In Fig. [3l we have shown how the two switch model proposed in [|5) can 
be expressed as a special case of the model we have proposed. There are three possibilities, 
viz (i) the CR transmitter and the receiver detect the same channels to be free. Here the CR 
transmitts using perfect knowledge and this is the case of ideal detection i.e. Pf a — Pmd = 0, 

(ii) the CR tranmsitter finds certain channels to be free which are not found free at the CR 
receiver hence if the CR transmits then it interfers with the PU which implies pf a ^ 0, and 

(iii) the CR receiver determines free channels but the CR transmitter does not which leads to 
a loss in transmission opportunity an hence, p m d ^ 0. 

III. Rate Region and Spectral Efficiency Factor for Ideal Detection 

In this section we present the rate region [0 and spectral efficiency factor for ideal detection. 
These results help in presenting the non ideal case in proper perspective. We then define 
spectral efficiency factor as the ratio of the sum capacity of the CR and PU to the capacity 
of PU in absence of CR. We discuss its variation as a function of p for different values of 
relative transmit power levels of CR and PU with help of an example. 

The system model for the ideal case is described in subsection III-AI The average capacity 
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for the PU, C p , and CR, C c , is given by flU 



C P = E M log^l + 
C c = E w [log ( 1 + 




) 



) 



bits/complex dimension. 



bits/complex dimension 



(4) 



The corresponding rate region is given by 



C = {(R c , R p ) |0 < R c < P C C , < Rp < (l-p)C p , 0<p< 1} 



(5) 



where the rates R c of the PU transceiver pair and R p of the CR transceiver pair is achieved 
through ideal white space filling [3]. Fig. |4] shows the rate region. 

The sum capacity of PU and CR is pC c + (1 — p)C p . There are three possible cases viz. 
(i) C p > C c , in this case p = maximizes the sum capacity, (ii) C p < C c , in this case p = 1 
maximizes the sum capacity, and (iii) C p = C c , in this case any value of p would maximize 
the sum capacity. 

We now define spectral efficiency factor in terms of the sum capacity as 

Definition 2. Spectral efficiency factor, r], is defined as the ratio of sum capacity of CR and 
PU to the capacity of PU in the absence of CR. Formally, 



• Dp = PUs data rate achievable with CR present, 

• D^ R = CRs achievable data rate 

• D p = PUs data rate achievable without CR present. 

With this definition, we have 
Theorem 1. The Spectral efficiency factor for ideal detection with average power constraints 



as P c and P p for the CR and PU respectively is given by 

where C c and C p are as given in (0). 

Proof: This follows from definition |2] and ©. ■ 
The value of this factor is always greater than unity as the CR transmits when the primary 
is idle, making more use of resource. It is of interest to see the variation of r\ with p. 

1 C c 



Corollary 1. ry increases with p with a rate of 



(i-pfc; 



Proof: As p increases, then (1 — p) decreases which implies that increases. Hence 

1 — p 

p 

increases with p. Since C„ and C c are independent of p, rj increases with p. By 

1 — p 

differentiating 77 w.r.t. p, we get the slope of the curve of r\. ■ 
We look at an example where fading is Rayleigh to better understand the consequence of 
Corollary [T] 

Example 1. Consider, the fading coefficients h p and h c to be uncorrelated and i.i.d over 

time, circularly symmetric zero mean complex Gaussian random variables with unit variance, 

denoted as h p , h c ~ CA^(0, 1). Thus \h p \, \h c \ are Rayleigh distributed with parameter value 

1, denoted as \h p \, \h c \ ~ Rayleigh{l). We assume a fixed P p = 10000 mW and we take 4 

p 

different values of P c such that RS = 10 log = {0,10, 20, 30} dB. We make plot of the i] 

vs p for these 4 different cases. Note that throughout the paper RS has the same definition 

unless otherwise specified. 

We observe from Fig. \5\ that for higher values of p, r\ increases rapidly irrespective of RS. 

1 

RS does not effect p. Hence, even if p increases slightly 77 increases rapidly. Looking at Fig. 
13 for RS = 10 dB, when we move from p = 0.8 to p = 0.98, the spectral efficiency factor 
increases from 5 to around 50 which is a ten times increase. Comparing this to RS = 



This is clear from Corollary [7] since — — increases rapidly, faster than increase in p and 



to 

dB, the growth is more rapid from p = 0.8 to p = 0.98 there is 20 times increase. Hence, 
for higher RS the growth is more rapid. To study the effect of r] for small value of p, i.e. 
p G [0,0.5] we blow up a portion of Fig. 4 in Fig. 5. Observe when RS is less, irrespective 
of p, the growth rate of rj is much higher than the case when RS is higher. This is evident 
from Corollary [Jj since lower RS would imply higher —. 

Note that while r\ for the ideal case is greater than unity, in the non ideal case it may not 
always be greater than unity as we see in the Section [V] 

Note that all the Figs. l8l - ITOl and [T3ll22l are plotted for Rayleigh fading channels with unit 
second moment. That is all h c , h p , h cp and h pc are complex Gaussian random variables with 
zero mean and unit variance. However, since we consider average or ergodic capacity, the 
analysis and conclusions are valid for any other distribution for which expectation is finite. 
Also note that the average received SNR at the PU receiver is assumed to be OdB unless 
otherwise specifically mentioned. 

IV. Rate Region and maximum sum capacity for Non Ideal Detection 

In this section, we find out the rate region for the non ideal case for which the system 
model is discussed in section III-Bl By considering all the cases defined in the time sharing 
random variable Q in definition [TJ we bring out the capacity expressions and then specify the 
rate region. We look at the effects of pf a , p md and (pf a ,Pmd) together on the rate region with 
the help of plots. We see that the rate region approaches that of the ideal case as (pf a ,Pmd) 
tend to zero. We then write the sum capacity and maximize it under two specific type of 
constraints viz, when p is specified and the maximization is over (pf a ,Pmd)- The other case 
is when {p fa-, Pmd) are specified and maximization is over p. 

We have, 

Lemma 1. The average capacity C and C' c of the PU and CR for average transmit power 
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constraints P p and P c are respectively as follows 



(l-p)((l-p fa )A p + Pfa B p )) 



(6) 



C' c 



P(l - Pmd)A c + (1 - p)p fa B, 



(7) 



where, 




The units are again bits/complex dimensions. 

Proof: We consider the 4 different situations described by Q. 
• For Q = 1 we have the baseband equations as given in © and reproduced here, 

Y c = h c X c + Z c . 

Hence, the capacity of CR given Q — 1 is 



Note that C p |(q = i) = 0. 
• For Q = 2 we have the baseband equations as given in © and reproduced here, 



Here, both the PU and CR treat each others transmissions as Gaussian noise with mean 
and variance E[|/i pc | 2 P p ] = P p and E[|/i cp | 2 P c ] = P c respectively. Hence, the total noise 
variance for the PU and CR is P p + a 2 and P c + a 2 respectively. Hence, capacity of the 
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PU and CR is 



Cp\(Q=2) = E| hj> | 



h \ 2 P 
log I 1 1 1 pl p 



> C C \(Q=2) = E|, 



log 1 + 



\h c \ 2 P c 
P P + <J 2 



For Q = 3. We have C c | (<9=3) = C p \ iQ=3 ) = 0. 

For Q = 4 we have the baseband equations as given in equation © 

hpXp -\- Zp. 

Hence, the capacity of PU given Q = 4 is 



a 



p\(Q=4) - E|h p 



\h \ 2 P 



Note that C c |(q =4 ) = 0. 
Finally we have, the capacity for CR and PU as 



c c - y]c c \( Q=i )P(Q 



i=l 



= (1 -p)((l-p /a )^+p /OJ B p )) 

4 

c; =^c p | (Q=i) p(g = z) 



i=l 



P(l - Pmd)-4 C + (1 - P)PfaB c 



(8) 



(9) 



These follow from the distribution of the time sharing random variable defined in the system 
model in subsection HI-Bl ■ 

Note that the values of A p , B p , A c and B c as specified in the Lemma [Q will be used 
throughout the paper unless otherwise specified. 

The corresponding rate region is given by 



C(Pfa,Pmd) = { (Rc,Rp)\0 <R P < (l-p)C 2 ( Pfa ), 

< Rc < (1 - p) C[ (Pfa)+ pCi {p md ) , < P < 1} 



(10) 
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where the rates, R p of the PU transceiver pair and R c of the CR transceiver pair are achieved 
through non ideal white space filling. Also, Ci{p md ) = (1 - p md )A c , C[(p fa ) = p fa B c and 
C 2 (pfa) = (1 - p fa )A p + p fa B p . Note that, 

C(Pfa,Pmd)\(p fa = Pmd =0) = C. 

where C is the rate region of the ideal case, specified in [5] In other words, C — > C as Pf a ,p md 
tends to zero. This implies that for smaller values of (pf a ,Pmd), the rate region in the non 
ideal case comes closer to that of the ideal case. Consider the rate region shown in Fig. [7] 
the interior of the red diagonal line along with the axes, specifies the rate region of the ideal 
detection and interior of the blue line along with the axes, specifies the rate region for the non 
ideal detection. Clearly, the non ideal case is the sub-region of the ideal case as it completely 
lies inside the ideal case region. Moreover, the cut on the vertical (or PU) axis, given by 
C 2 (pf a ), is purely due to pf a , while that on the horizontal (or CR) axis is due to p md and pj a 
both (given by C\{p md ), C'^pfa) respectively) as justified from the expression of C and C' c 
specified in Lemma [TJ Also note that for a R c < C[(pf a ) there is no decrease in R p . 

Fig. [8] shows that the rate region for different values of pf a when the p m d is constant. Note 
that the non-ideal rate region approaches the ideal rate region as pf a decreases, along the 
Y axis only. This is because the pf a effects the PU more than the CR. Fig. [9] shows that 
the rate region for different values of p md when the pj a is constant. Observe the change in 
rate region is only on the X axis. Fig. [TO] gives the rate region when both pf a and p md vary 
simultaneously. 

A. Maximum sum capacity 

The sum capacity of CR and PU is given by 

C' c + C' p = p(l - p md )A c + (1 - p) Pfa B c + (1 - p)((l - Pfa )A p + Pfa B p ) 

= -A c pp md -(l-p)(A p -B p -B c ) Pfa + pA c + (l-p)A p . (11) 
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The maximum sum rate that the CR and PU can transmit, for a given p, can be found out 
by solving the following optimization problem. The data rate achievable by the CR is 



p.f 



max <^ - A c pp md - (1 - p)(A p - B p - B c )p fa + pA c + (1 - p) A p \ , 

fa,Pmd { ) 

s.t, , < p fa < 1, < p md < 1. (12) 

It is easy to verify that the maximum is attained when p/ a = Pmd = 0. The value of sum 
capacity will then be equal to that in the ideal case. 

The inverse problem of the maximum sum rate is when we have a given pj a and p m d and 
we want to see what value of p maximizes the sum capacity. 



max | - A c pp md - (1 - p)(A p - B p - B c )p fa + pA c + (1 - p)A p j, 

s.t,0<p<l. (13) 

The objective function can be re-framed as 

max j((A p - B p - B c )p fa - A c p md + A c - A p )p - (A p - B p - B c )p fa + A p |, 

si,0<p<l. (14) 

The solution to this depends upon the term M = {(A p — B p — B c )pf a — A c p md + A c — A p } 
and is given by 

1, if M > 

P= I 0, if M < • (!5) 

any value in [0, 1] if M = 

V. Spectral Efficiency Factor in the Non Ideal Detection 

In this section, we use expression of the sum capacity of CR and PU to express spectral 
efficiency factor in terms of p, p md and pj a . We look for the pairs (pf a ,p m d) which result in 
spectral efficiency factor being more than unity for a given value of p corresponding to the 
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case where the sum capacity in the CR system is at least as much as was the capacity of PU 
alone without CR. Then those (pf a ,Pmd) pairs are explored which guarantee that the PU's 
performance loss is within a limit. All such pairs make up admissible regions. Note that a pair 
(pf a ,Pmd) that is not admissible is termed as inadmissible and all such pairs constitute the 
inadmissible region. Appropriately, we define the two types of regions (i) where the rj > 1 - 
called weakly admissible region and (ii) where the PU does not loose beyond a limit - called 
strongly admissible region within a specified loss factor for PU. We also find what pairs that 
are admissible if the PU cannot incur any loss and term these as strong admissible pairs. We 
then discuss the effects of pf a , p md and p and show that pj a has more profound effect than 
p m d on the performance of CR. 

The spectral efficiency factor for the non ideal case is as shown below. 

Theorem 2. The spectral efficiency factor in the non ideal case % with average power 
constraints P p and P c for the PU and CR respectively is given by 



where C' c and C are specified in Lemma U\ 

Proof: Looking at definition [2] and the definition of the achievable data rates in the non 
ideal case from Lemma \T\ we have the result. ■ 
In the ideal case the spectral efficiency factor depended only on one parameter p, while in 
the non ideal case here it depends upon pf a and p m d as well. Hence the following corollary 
helps in studying the rate at which spectral efficiency increases/decreases with change in any 
one of these parameters keeping the others fixed. 

Corollary 2. Spectral efficiency factor n increases with p for a fixed p m d with a rate given 



by 




(17) 
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Proof: By partially differentiating (TT6T ) w.r.t p we get the result. ■ 

Corollary 3. Spectral efficiency factor n decreases with p md for a fixed p with a rate given 
by 

pA c 



(i-p)V 

Proof: By partially differentiating ([Tol w.r.t p md we get the result. 

Corollary 4. Spectral efficiency factor ft decreases with pf a with a rate given by 

(Ap — B p — B c ) 



(18) 



(19) 



Proof: By partially differentiating (TT6T ) w.r.t pf a we get the result. ■ 
It is interesting to note that the rate of increase in f) with p is independent of pf a . Also the 
rate of decrease in f) with p md is independent of pf a . Moreover, the rate of decrease in fj with 
Pfa is independent of both pf a and p md . Hence, we can say that pf a decreases the spectral 
efficiency factor of the system equally at all occupancies. 

Observe from Theorem |2] and Lemma [TJ that not all pairs of (pf a ,Pmd), will result in the 
spectral efficiency factor to be greater than one. This is because if there is high pf a then due 
to interference both CR and PU will have reduced transmission rates. If there is high p md 
then the CR will have less transmission rate. To study precisely the effects of (pf a ,Pmd) and 
to characterize them in a region for which f) > 1 we define the notion of weak admissibility 
as follows 

Definition 3. A pair (pf a ,Pmd) is said to be weakly admissible for a Bernoulli occupancy 
p and average power constraints P p and P c on the PU and CR respectively, if the spectral 
efficiency factor f) > 1. 

When the CR interferes in the PU communication, there will be a loss to the PU. Suppose 
we are given a limit beyond which the PU cannot incur a loss. It is of interest to see what 
values of (pf a ,Pmd) can ensure that the loss to the PU is below this limit. We model this in 
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the form of the loss factor. The loss factor of the PU is the fraction of data rate it looses due 
to the intervention of the CR in its communication. Formally, 

Definition 4. The loss factor of the PU is 

g 

Note that < 7 < 1. We now wish to see what (pf a ,Pmd) will guarantee a loss below 7. 
This gives rise to notion of strong admissibility with a loss factor 7 as follows, 

Definition 5. A pair (pf a ,Pmd) z * sa id to be strongly admissible with a loss factor 7 for 
a Bernoulli occupancy p and average power constraints P p and P c on the PU and CR 
respectively, if C > 7(1 — p)C p . In particular if 7 = 1 we say that (pf a ,Pmd) are strongly 
admissible. 

Now, we characterize the admissible regions based on the above definitions. 

A. Characterization of Admissible Regions 

In this subsection we characterize three types of admissible regions viz (i) weakly admis- 
sible region - Theorem |3] (ii) strongly admissible region with a loss factor 7 - Theorem @] 
and (iii) strongly admissible region - Theorem |5] 

Theorem 3. The pair (pf a ,Pmd) are weakly admissible if 

(1) 0< Pfa < 1. 

n\ ^1 f l ~P\ ( Ap~B p -B c \ 

(2) p^Kl-l—j { j p fa . 

Proof: The first constraint is kept there to realize that the other constraint can be satisfied 
for almost all values of (pf a , Pmd) if Pfa is not confined to the interval of [0,1]. From definition 
of weak admissibility 

fj > 1 or, 

p 1 Pfa 

B p )) + p(l - p md )A c + (1 - p)p fa B c > (1 - p)A p (20) 
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Grouping terms of pf a and p m d together, we have 



(1 - p) [A p - B p - B c ]p fa + pA c p md < pA c or, 



(l-p\ ( A p -B p -B c \ 



Observe that the region for strong admissibility with a factor 7, only depends on pf a . This 
is because strong admissibility is concerned with the data rate of PU only. For a pf a that is 
strongly admissible with a factor 7 all values of p m d are admissible. Therefore we plot pf a 
for strong admissibility. Also we provide characterization of strong admissible regions with 
a loss factor 7 in terms of pf a only as follows 

Theorem 4. The value of pf a is strongly admissible for a given loss factor 7 of the PU if, 




Proof: We have, 7 < 1. From (1), Lemma [Hand definition @] we have 



(1 -p fa )A p +p fa Bp < jA p or, 




(21) 



Now since, pf a < 1, the result follows. 



Theorem 5. Strongly admissible region is characterized by pj a = 0. 



Proof: For strong admissibility we have 7=1. Thus, 



(1 - p)((l - p fa )Ap + p fa B p )) = (1 - p)C } 



p 



(22) 
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Note that, C p = A p , hence we have 

p fa (A p - B p ) = 0. (23) 

A p = B p in the limit would imply that P s <C a 2 which is never the case. Hence, pf a — 0. ■ 
Theorem [5] implies that any detection technique with pj a > will lead to spectral efficiency 

loss for the PU. Since all the practical techniques, will have pj a > 0, we have introduced the 

notion of strong admissibility with a loss factor. 

Note that the weak admissible region occupies a non zero area in the (pf a ,Pmd) plane for 

p > as evident from Theorem [3] An instance of the weak and strong admissible region is 

shown in Fig. [TTJ The set of strongly admissible (pf a ,Pmd) pairs (i.e. 7 = 1) is a sub-region 

of the weakly admissible pairs (pf a ,Pmd) as shown in Fig. [TTJ 

B. Strong Admissible pj a with a loss factor 7. 

A typical curve of the strongly admissible pf a with a loss factor 7 against 1 — 7 is shown 



in Fig. [T2J As we decrease 7, 1 — 7 increases the rate at which the admissible value of pj a 

increases is ^ p — . 

Ap — Bp 

Theorem 6. The range of strongly admissible values ofpf a for a loss factor 7 is independent 

B„ 

of the PU channel occupancy. Furthermore, for a loss factor of 7 < — the range of values 

Ap 

is the entire interval [0,1]. 

Proof: From Theorem [4]it is clear that admissible pj a values for a given 7 is not depen- 
dent on p. The second statement of the Theorem, follows from the Theorem [4] substituting 

Pf a — 1 in the second inequality. ■ 

B v 

We call 7 = -p as the full admissible point. Observe that for 7 = 1, i.e. the strong 

A p 

admissible case p/ a = 0. Also note that the admissible values of pf a are highly dependent on 
the relative power levels of the PU and CR (RS) transmission. As the value RS increases 
the values of strongly admissible pf a for a given loss factor also increases as shown in Fig. 
Similarly as RS increases the value of the full admissible point increases as illustrated 
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in Fig. [141 

Next we consider the weak admissible region. Fig. [15] shows the spectral efficiency factor 
plot for a particular value of p and RS with the admissible and inadmissible regions marked. 
Fig. [ID shows the spectral efficiency factor for various values of p for a constant RS. Observe 
that the admissible region increase as p increases. This implies that for channel with low 
primary occupancy a lousy detector (with high Pf a ,Pmd) w iU a l so result in increase in spectral 
efficiency. Similarly for pf a = the PU will transmit its original capacity without CR and 
hence any value of p rnd would be weakly admissible. Also observe from Fig. [20] that when p 
is fixed and RS is increased, the admissible region increases. This is justified because, when 
the RS is more the CR transmits relatively lesser power which causes lesser interference 
when compared to when RS is lower. This translates into higher values of admissible pf a for 
higher RS when compared to the lower RS. 

We show the weakly admissible and inadmissible regions in the Figs. [17] and [T8l 

VI. Comparison of Various Detection Techniques in light of the Admissible 

Regions 

We compare few of the standard detection techniques viz. energy detection (ed), matched 
filter (mf) and magnitude squared coherence (msc). 

For every detection technique (dt) we have a relation fat between the p m d and pf a , i.e. 
Pmd = fdtipfa) for a given average received signal to noise ratio (SNR). The function f dt for 
all standard detection techniques is monotone non increasing. For a given pf a we wish to see 
whether any value of p m d is admissible. 

J l— — -J pf a . We then wish to find for a detection 

technique whether for a given pj a , fdt(Pf a ) < h(pj a ) holds. If yes, then we say that the 
detection technique operating at that (pf a ,Pmd) pair is admissible and consequently all such 
points on the curve of detection technique are referred to as admissible region of the detection 
technique. Note that, this admissible region of a detection technique is dependent on average 
received SNR. 
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We provide the functions f dt for three detection techniques mentioned above. In what 
follows the value of p m d is a function of pj a for all other parameters fixed. 

Pmded = fediPfaed) = Q x 2 (^P^ 1 ^1 - p faed , L, ^f P j V for a fixed Af,L and P p , 

Pmdmf = fmf(Pfamf) = 1 ~ Q \Q~ l {Pfamf) ~ ' for a fixed E and a ' 

Pmdmsc = fmsc(Pfamsc) = PcDF ( U ~ P/amslJ l L > ItN , f ° r a fixed L and 7- (24) 

where, Q x 2(x, v,8) is the non central x 2 distribution of x with w degrees of freedom and 
positive non centrality parameter 5 and P _1 (x, a) is the inverse lower incomplete T function. 
Q(x) is the standard Q-function with Q^ 1 being the inverse Q function. Pcdf(\^\ 2 \L, |7| 2 ) is 



as given in IIT9I . Here, L denotes the number of disjoint sequences that an original sequence 

N 

of length TV is divided into each of length M so that M = — . E is the energy of the 

Lj 

transmitted signal. I7I 2 is the magnitude squared coherence which is the magnitude squared 



of the spectral coherence. For more details, the reader is referred to [1201 and the references 
therein. 

We now plot p m d vs. pf a and look for those areas where the curve of detection region lies 
inside the weak admissible region. For all points lying on that part of the curve the detection 
technique is admissible. 

Figs. [2T1 [22] show that, for low SNR regimes, all values of (pf a ,Pmd) offered by the matched 
filter detector are admissible. For the magnitude squared coherence detector majority of the 
points on the curve are admissible except those at higher end of p/ a values. For the energy 
detector many points on the curve lie outside the admissible region. Hence, using energy 
detector, one has to be careful about the operating point (pf a ,Pmd)- It follows that the matched 
filter performs better than magnitude squared coherence detector which performs better than 
the energy detector. These observations are in line with those in [BTTl . 
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VII. Conclusion 

In this paper we have capacity region of a point to point CR channel. For ideal detection 
of PU occupancy we have seen the limits on rates. We have defined spectral efficiency factor 
and studied its variations with respect to the occupancy probability for different relative power 
levels of the PU and CR. Then we built the case where the CR performs non ideal detection 
of the PU presence and specified capacity region for the same. With expressions derived we 
brought out the limits, under different occupancies and relative PU and CR power levels, on 
the false alarm and missed detection for which the spectral efficiency of the overall system 
increases and also for which the PU is not at a loss greater than a specified value. We have 
discussed the effects of false alarm, missed detection and channel occupancy with respect to 
spectral efficiency factor. Finally we compared our result with standard detection techniques 
viz. energy detection, matched filter and magnitude squared coherence and found that the 
matched filter performs best followed by magnitude squared coherence followed by energy 
detection with respect to spectral efficiency factor. 
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Fig. 1. The Cognitive radio channel when CR ideally detects the PUs presence. 




Fig. 2. The Cognitive radio channel when CR non ideally detects ideally detects the PUs presence. The dashed instance 
shows when PU transmits alone and dotted when PU and CR transmit together and solid when CR alone transmit. 
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Spectral hole detected at CR Txt and Rx both the case of 
single switch model, p/ a — Pmd = 




Spectral hole detected at Spectral hole detected at 

CR Txt but not at CR Rx, CR Rx but not at CR Txt, 

Pfa 7^ Pmd ± 



Fig. 3. The two switch model shown as the special case of the non ideal case system model. 




Fig. 4. Rate region of ideal CR. 
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Fig. 5. Dependence of r\ on p for various relative power levels of PU and CR. 





Fig. 8. Rate region of non ideal case compared with non ideal for fixed p m d but variable pf a . 
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Fig. 9. Rate region of non ideal case compared with non ideal for fixed pf a but variable p m d- 




Fig. 10. Rate region of non ideal case compared with non ideal for two pf a and p mt j values. When pf a and p m d both 
decrease the rate region becomes closer to the ideal case. 
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Strong admissibility line pj a = 




A>(l-7o) (1,0) 
— Bp 

Fig. 11. Weak admissible region appears as shown above. It increases with increase in value of p and when PU to CR 
power ratio increases. Note here the value of pf a which causes the change in admissible region. This is explained later. 
Strong admissible region for all values of p and CR and PU powers transmit is the line pf a — 0. Strong admissible region 
with a loss factor 70 is the shaded rectangle. As 70 increases the size of the rectangle decreases and eventually becomes 
the line pf a = for 70 = 1. Note that strong admissible with a loss factor 70 is not necessarily a sub region of the weak 
admissible region. 




Slope 



A p — Bp 



1-7 



Fig. 12. pfa vs (1 — 7). As (1 — 7) increases, the maximum value of pf a increases up until 7 = — — , where pf a < 1 
becomes admissible. 




Fig. 14. Plot of full admissible point vs relative power levels of the CR and PU; RS. for p = 0.4 




Fig. 16. Spectral efficiency factor variation vs pf a and p m d for various values of p and RS = dB, highlighting the 
change in admissible region with p. A cut has been introduced to separately depict the spectral efficiency factor in the 
weakly admissible regions only. 
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Admissible regions 



p 6 04 



Fig. 17. pf a vs p m d plot showing weakly admissible regions. Note that the weakly admissible region for p = po (say 0.6) 
is the sum of shaded weakly admissible regions for all p less than po. So the admissible region for p — 0.6 is the sum of 
shaded weakly admissible regions for p = 0.2 and 0.4 plus the region shaded for p = 0.6. 




Fig. 18. pfa vs Pmd plot showing inadmissible regions. 




Fig. 20. pf a vs pmd plot for a fixed p = 0.5 for various values of relative power levels of PU and CR, RS. 
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Energy detection SNR = -24 dB 

Matched filtering SNR = -24 dB 

Magnitude Squared coherence SNR = —24 dB 
Admissible regions SNR = —24dB 
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Fig. 21. Admissible regions for p = 0.2 and RS = 20dB for the energy detector, magnitude squared coherence and 
matched filter. The average received SNR is — 24dB 











Energy detection SNR = -26 dB 

Matched filtering SNR = -26 dB 

Magnitude Squared coherence SNR = -26 dB 
Admissible regions SNR = -26dB 
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Fig. 22. Admissible regions for p = 0.2 and RS = 20dB for the energy detector, magnitude squared coherence and 
matched filter. The average received SNR is —26dB 
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